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Qh! Abstract 

^ ■ 

Many extensions of the standard electroweak model Higgs sector suggest that 

the main Higgs decay channel is "invisible", for example, h — > J J where J denotes 

k> I the majoron, a weakly interacting pseudoscalar Goldstone boson associated to the 

j_j ■ spontaneous violation of lepton number. In many of these models the Higgs boson 

Cd ' may also be produced in association to a massive pseudoscalar boson (HA), in 

addition to the standard Bjorken mechanism (HZ). We describe a general strategy 

to determine limits from LEP data on the masses and couplings of such Higgs 

bosons, using the existing data on acoplanar dijet events as well as data on four and 

six b jet event topologies. For the sake of illustration, we present constraints that 

can be obtained for the ALEPH data. 
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1 Introduction 

The Higgs particle remains one of the missing hnks in the otherwise well tested standard 
electroweak model (SM) [l|]. Its mass is not fixed theoretically and experimental efforts 
to search for the higgs boson, notably at LEP, have produced a lower bound ~ 63 GeV 
on its mass. This bound applies to the standard model (SM) Higgs. 

There exist well motivated extensions of the SM which are characterized by a more 
complex Higgs structure than that of the SM P, ^. These include (i) the minimal super- 
symmetric standard model (MSSM) {ii) generic two or more Higgs doublet models and 
{hi) the majoron type models characterized by a spontaneously broken global symmetry. 
This symmetry could either be a lepton number [^ 0, a combination of family lepton 
numbers or an R symmetry in the context of supersymmetry [0. These extensions con- 
tain one or more parameters in addition to the Higgs mass. Nevertheless, a considerable 
region of the parameter space has already been ruled out in these models using the LEP 
data §. 

Among the three extensions mentioned above, the majoron type models are quali- 
tatively different compared to the other two as well as to the SM. These models contain 
a massless goldstone boson, called majoron. There are many types of majoron models 
0]. Many of such models are characterized by the spontaneous violation of a global U{1) 
lepton number symmetry close to the electroweak scale. This could also have important 
implications for the structure of the electroweak phase transition and the generation of 
the electroweak baryon asymmetry [^. 

In such models the normal doublet Higgs is expected to have sizeable invisible decay 
modes to the majoron, due to the strong higgs majoron coupling []10[ - |T^. This can have 



a significant effect on the Higgs phenomenology at LEP. In particular, the invisible decay 
could contribute to the standard signal looked for at LEP namely two acoplanar jets and 
missing momentum. This feature of the majoron model allows one to strongly constrain 



the Higgs mass [^, |T6|, |T^ in spite of the occurrence of extra parameters compared to 
the standard model. In particular, the limit on the predominantly doublet Higgs mass 
can be seen to be close to the SM limit irrespective of the decay mode of the Higgs boson 

n, nn, mi- 



Apart from giving an interesting twist to the Higgs search strategy, the majoron 



models are quite well motivated theoretically 0. They allow the possibility of spon- 
taneously generating lepton number violation and hence neutrino masses. The Higgs 
couplings to majorons in these models were discussed at length in |rS|. 



As far as the neutral Higgs sector is concerned, all the models fall in two categories. 

(A) Minimally extended SM characterized by the addition of a singlet Higgs: Such models 

are typical examples of mechanisms where the neutrinos obtain their mass at the tree level 

T9| . (B) Models containing two Higgs doublets and a singlet under the SU{2) x f/(l) 



group: These models though more complex arise naturally in trying to generate the 
neutrino mass radiatively. The Zee type |^ model modified to obtain spontaneous lepton 
number violation is a typical example |^ in this category. These models have the virtue of 
explaining smallness of neutrino masses without invoking any high scale. Moreover, type B 
models are different from the type A ones as far as the Higgs phenomenology is concerned 
since they contain a massive pseudoscalar which can be produced in association with the 



Higgs as in MSSM. All the existing analysis of the invisible Higgs search [JS], |T^, |17| have 
concentrated on type A models and hence on the Bjorken process. 

In this letter we extend this analysis to include the type (B) models and particularly 
the associated production of the pseudoscalar Higgs boson. We describe a general strategy 
to determine limits from LEP data on the H and A masses and couplings using the existing 
data on acoplanar dijet events. We also show how one can set decay mode-independent 
limits on H and A masses and couplings by including data on four and six b jet events. For 
the sake of illustration, we display the constraints that can be obtained for the published 
ALEPH data. 



2 Higgs Boson Production and Decay 

Let us first recall the salient features of the type (B) models. The Higgs sector is charac- 
terized by two doublets 01^2 and a singlet a. The part of the scalar potential containing 
the neutral Higgs fields is given in this case by 

Ai2(0l0l)(0|02) + Ai3(0l0l)(aV) + A23(0^02)(^V) 

+5(4^,)(4^,) + ^n[i4<p2f+ h.c] (1) 



where a sum over repeated indices i=l,2 is assumed. Here 0i^2 are the doublet fields and 
o corresponds to the singlet carrying nonzero lepton number. 

In writing down the above equation, we have imposed a discrete symmetry 02 -^ —<p2 
needed to obtain natural flavour conservation in the presence of more than one Higgs 
doublet. For simplicity, we assume all the couplings and VEVS to be real. 

After minimization of the above potential one can work out the mass matrix for the 
Higgs fields. To this end we shift the fields as (i=l,2) 

Ri + Hi 



Vi 



a 



v^ V2 



V3 R3 + ih 



V2 V2 

The masses of the CP even fields Ra (a=1...3) are obtained from 
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\viVYVz A23t^2^^3 2A3t;| 

The physical mass eigenstates Ha are related to the corresponding weak eigenstates as 



Ha — Oab Rb 



where, O is a 3x3 matrix diagonalizing Mj^ 



O M^O^ = diag (M^, M^, M| 



(6) 



(7) 



It is convenient to parametrize the matrix O as 

/ 



O 



^0.^0 ^a^uj ^a^uj^O ^a^uj Ca^uj^O 



^a^9 ^ofiijj ^OL^LO ^6 ^a ^lo ^a ^lo ^6 



se 



SujCg 



CujCe 



where c„ = cos a etc. 



In addition, there exists also a massive CP odd state A, related to the doublet fields 
as follows 

A=^iv2h-vj2). (9) 

Its mass is given by 

Ml = -kV'^ (10) 

where V = {vf + vlY^'^. Thus, after spontaneous SU{2) x f/(l) x U{1)l breaking, we have 
a total of three massive CP even scalars H^ (i=l,2,3), plus a massive pseudoscalar A and 
the massless majoron J, simply given as J = I3. 

At an e~^e~ collider there are two main production mechanisms that allow the pro- 
duction of the higgs boson through their couplings to Z. The relevant couplings for their 
production through Bjorken process (HZ) are given as follows (a=1...3) 

Chzz = {V2GFf'^Ml Z^Z^i'^O.a + '^02a]Ha. (H) 

As long as the mixing appearing eq. (^TJ) is O (1), the Higgs bosons can have significant 
couplings and hence appreciable production rates through the Bjorken process. 

In addition, the Ha can also be produced in association with the CP odd field A 
through the HA coupling 



^HAZ — 7, — Z^ 

cosdw 



yOra - y02a^ 



Ha d^ A. (12) 



In what follows we assume that at LEP only the the lightest of the CP even Higgs 
boson, denoted hj Hi = H can be accessible. Using the matrix O of eq. (H) in eq. ([TI|) 
and eq. (|12D one gets the couplings of H to ZZ and HA as 



Chzz = {V2GFf''^Ml Zf.Z^'cosiP - a) cosOH (13) 

where tan /3 = — and 

Chaz = ^Z^'sm(l3-a)cose H d^" A. (14) 

cosOw 

As expected, one recovers the well known expressions for the two doublets case 
in the limit 6^0. In particular, note that for sizeable value of cos 6', one cannot 



simultaneously suppress the production of Higgs through eq. ([T3|) and eq. (|1^) if it 
is allowed kinematically. 



Now we turn to Higgs boson decay. For Higgs boson masses accessible at LEP 
energies the main CP even Higgs decay modes are into bb and JJ. The couphng of the 
physical Higgses to J follows from eq. (|I]) . One can express this coupling entirely in terms 
of the masses M^ and the mixing angles in the matrix O 

Cj = -J\2\3V3R3 + X13V1R1 + X23V2V3R2) (15) 

= ^^{Ml)3aRa (16) 

= ]^{V2GFf'H^n^{0^)3aMlHJ^ (17) 

where tan7 = — . We have made use of eq. (^ and eq. (|^) in writing the last line. 



Note from eq. (p!5| ) that the CP even scalar Higgs bosons Hi couple strongly to a pair 
of majorons leading to the invisible decay signature. In contrast to ifj, the pseudoscalar 
A does not decay into one or three majorons since the couplings AJ^ or AH J, although 
possible in general, do not exist at tree level in our simplest model described above. 
Because of the form of the scalar potential given in eq. (|I|) only even powers of the majoron 
field J = R3 appear once the expansions in eq. (0) and eq. (^ are used. As a result 
the couplings AJ^ and AH J are absent from the potential and therefore the A decays 
visibly to fermion antifermion pair 0. The branching fraction B^ for A —* bb is nearly 
one. Deviation from unity of this branching ratio is model dependent. If all fermions 
obtain their mass by coupling to only one higgs boson, as in the models considered in ref. 
jTSl , H, then one has 

B. = ^ (18) 

where 

mjjl - Amj/miy/^ 

""^ y ml{l-Aml/m\y/^ ^ ' 

and the sum is over all the fermions except b. The rate for H ^ bb can be expressed as 

T{H -^ bb) = ^^^^ Mgmg(l - Aml/M^f^ cos^ a cos^ 9 (20) 

On the other hand the width for the invisible H decay can be parametrized by 

T(H -. JJ) = y^Mf, (tan7 On)' (21) 



**Note that A could decay invisibly at the tree level through the process A -^ Z*H -^ vvJJ. However 
the branching ratio in this mode is quite small. 



For many choices of parameters the invisible decay mode can be rather important and, in 
fact, provide the strongest hmits. Finally, note that additional decay modes of H to A A 
may also exist. 

3 Search Strategy and Analysis 

Since A can decay only visibly, one expects dijets + missing momentum as a signal of 
the Higgs production in our simplest model. This signal arises from three processes - 
i.e. Z* -^ vv with E -^ bb, or Z* -^ qq with H^JJoiH^JJ with A -^ bb. For 
each process one has a sizeable missing momentum which is aligned neither along the 
beam nor along the jets. In contrast, for the SM background, the missing momentum 
arises from i) jet fluctuation (including decay u from b, c and r jets) in which case it 
is aligned along the jets, and ii) initial state radiation (ISR) or the two-photon process 
e+e^ -^ (e+e^)77 in which case it is aligned along the beam direction. This enables 
one to eliminate the SM background in the dijet + missing momentum channel by a 
suitable combination of kinematic cuts without depleting the signals seriously. In fact, 
this procedure has been extensively used to search for the SM Higgs signal (first process) 
in the LEP data and obtain the corresponding mass limit 0. More recently the analysis 
has been extended to the invisibly decaying Higgs signal (second process) and obtain the 
corresponding mass limits for the A type majoron models |l^, [l^, [l^. In this section 
we shall extend the analysis further to include the Higgs signal from the third process of 
associated production (HA) and study the resulting mass limits for the B type models. 

We shall use a parton level Monte Carlo event generator, which has been shown [|1^ 
to reproduce the signals obtained with the full Monte Carlo program ^ quite well. 

For our illustrative purposes, below we will determine the corresponding Higgs boson 
mass limits which can be obtained from the ALEPH data sample of ref. f^ based on a 
statistics of ~ 1.23 million hadronic Z events. A detailed account of the experimental cuts 
can be found in . We shall only summarize the main features. One starts with a visible 
mass cut, Mjj < 70 GeV, to ensure a sizeable missing energy (^) and momentum (p). A 
low angle cut, requiring the energy deposit within 12° of the beam axis to be < 3 GeV 
and that beyond 30° to be > 60% of the visible energy, removes jets close to be beam 
pipe where measurement errors can simulate a large ^. An acoUinearity cut, requiring the 
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angle between the two jets to be < 165°, suppresses the Z —>■ qq and t^t~ background 
where the ^ can come from jet fluctuation (including escaping u). Moreover an isolation 
cut on ^removes the ^ background from the fluctuation of any one of the jets. A cut on 
the angle (a) of ^ with respect to the beam axis, tana > 0.4, suppresses the background 
from ISR and two photon processes. The I^ background coming from the jet fluctuations 
along with an ISR are removed by an acoplanarity cut, requiring the azimuthal angle 
between the two jets to be < 175°. An acoplanarity cut for 3-jet like events, requiring 
the sum of the 3 dijet angles to be < 350°, removes the ^ background arising from the 
fluctuation of these jets. The remaining few events are residual two-photon events, which 
are removed by a total p^ cut. The cuts remove all the events from the data sample P], 
while retaining ^ 50% of the signal for each of the three processes mentioned above. 

We denote the number of signal events for the three processes, after the cuts, by 
NsM,Njj and Na respectively, assuming no suppression due to the mixing angles or 
branching fractions in each case. Then the expected number of signal events, after incor- 
porating these effects, is given by 

Ne,pt = e| [BNj + (1 - B)Nsm] + e\BABNA (22) 

where 

B = BR{H -^ J J) (23) 

is the branching fraction for the H decaying into the invisible mode and e^ = sin(/3 — 
a) cos 9 and e^ = cos(/3 — a) cos 6. 

The number of signal events for the three processes that pass the cuts are shown in 
Fig. 1. As we can see Na ^ Nsm, Njj thus implying that, if kinematically open and 
not suppressed by mixing angles, associated production tends to give the strongest limits. 
We shall now consider the limits which can be obtained by comparing eq. (22) with the 
95% CL limit of 3 events, corresponding to events in the data sample |§] after the cuts. 

First, we assume that B=0. In this case we get the limits corresponding to the SM 
decay mode. Note that this is the weakest possible limit for the Bjorken coupling strength 
e\. This is in complete agreement with the results of ref. |[T5| , |16[ as seen in Fig. 2. Note 
that in this case we can set no limits on e^. 

We can also obtain limits assuming B=l. From Fig. 1 one can see that, except for the 
region very close to the edge of the phase space, uia + t^h ~ '^z, the major contribution 



comes from the associated production [Na) allowing the corresponding coupling strength 
parameter e\ to be strongly constrained, as seen in Fig. 3. It is easily seen that the limits 
close to the edge of the phase space are typically at the level e\ ~ 0.05. However, very 
close to the kinematical limit the Bjorken process becomes important and could further 
strengthen the limits we have obtained. Outside the region where associated production 
is kinematically possible, the Bjorken process can still take place and the limits on e^ are 
the same as obtained in ref. [T3, |TB|, 0. The limits on the coupling strengths e\ and e% 



can be improved further by using more recent data from LEP. 

In the case B=0 considered earlier one does not get any limit on e\ since in this 
case the associated production does not lead to the dijets + missing momentum signal, 
we have considered so far. However one could derive limits on the e\ using data on 
four and six b jets topologies, assuming good b identification. As an illustration of what 
can be achieved we assume some typical value for the corresponding branching ratios 
BR{e~^e' ^ H A —>■ 46) < 5 x 10~^. Using this value we can determine the excluded 
region for different values of e\ as a function of rriA and mn from 

BRha = \ {BR),,\' el , (24) 

where BR^y is the Z decay branching ratio into one generation of neutrinos in the SM, 
the factor 1/2 refers to the scalar HA decay mode and A is the corresponding phase space 
factor 



[s + m?fj — m\y — As m\j 

X{s,mH,mA) = 

s 

The resulting excluded regions are shown in Fig. 4. On the other hand if H has a mass 
bigger than 2mA then eq. ( P^ would give a weaker limit since it should now be multiplied 
by the branching ratio for H into 2b. But now additional information can be obtained 
from the study of events with six b jet topologies. One may obtain a similar limit using 
an illustrative reference limit for the corresponding branching BR[e~^e^ -^ H A ^ 66) < 
4 X 10^^ 



4 Discussion 

Here we have shown that the production of invisibly decaying Higgs bosons in Z decays 
through the associated channel HA can lead to substantial limits on Higgs boson masses 
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and couplings. Such invisibly decaying Higgs bosons arise in a wide class of SU{2) ® U{1) 
niajoron-type models, such as the majoron extension of the Zee model for radiatively 
induced neutrino masses, as well as extensions of the minimal supersymmetric standard 
model with spontaneously broken R parity 0. 

Note that it is also possible to derive a truly decay mode independent limit on the 
H mass in the case of associated HA production by allowing B to vary from to 1 and 
by combining the data on dijet + missing momentum, with those from four and six b jet 
searches described above. 

It is important to update the limits obtained here by using the full statistics provided 
by the four LEP experiments and take them into account in designing the strategies to 



search for the Higgs boson at higher energies, such as at the LHC [£2[ and NLC ||2T |. 



Finally we note that a similar analysis may be performed in the context of the 
minimal supersymmetric standard model, where the invisible decay of the Higgs bosons 
into neutralinos may take place. A major difference in the latter case arises from the fact 
that in the MSSM case the pseudoscalar boson A is the one most likely ||2^ to decay 
invisibly, whereas in the case discussed in the present paper the pseudoscalar A always 
decays into bb. 
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Figure Captions 
Figure 1 

Numbers of expected dijet+ missing momentum events for the processes e~^e^ -^ Z*H -^ 
vv hi (NcjAf), e+e- -^ Z*H -^ qq J J (Njj) and e+e" -^ H A ^ J Jbb (N^a) after 
imposing ALEPH cuts [|, Q]. 
Figure 2 

Illustrative limits on e^ for visibly decaying H (B=0) from the Z*H — > z/z/ + two jets 
channel as a function of rriH using the ALEPH data of ref . [0 . 
Figure 3 

Limits on e\ in the mAfnH plane, based on the e"'"e^ ^ H A ^ J Jbb production channel. 
We have assumed BR {H —>■ J J) = 100% and a visibly decaying A with branching ratio 
into bb given in eq. (18) 
Figure 4 

Limits on e\ for the B=0 case based only on e+e^ —>■ H A ^ bbbb channel. We have as- 
sumed a hypothetical sensitivity for the 46 channel of BR(e"''e~ -^ H A —* 46) < 5 x 10~^. 
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